Dynamically Driven Renormalization Group Applied to Sandpile Models 
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The general framework for the renormalization group analysis of self-organized critical sandpile 
models is formulated. The usual real space renormalization scheme for lattice models when applied 
to nonequilibrium dynamical models must be supplemented by feedback relations coming from the 
stationarity conditions. On the basis of these ideas the Dynamically Driven Renormalization Group 
is applied to describe the boundary and bulk critical behavior of sandpile models. A detailed de- 
scription of the branching nature of sandpile avalanches is given in terms of the generating functions 
of the underlying branching process. 

PACS number(s):64.60.Lx, 05.40. +j, 05.70.Ln 



CO 



I 

C 

O 

o 



> 
in 
o 

(N 
(N 
O 
00 
Os 

C3 



-a 
c 

o 
o 



X 



I. INTRODUCTION 

In the last decade it has been recognized that a 
fair amount of physical phenomena are characterized by 
strong fluctuations and long-range correlation functions. 
According to the theory of equilibrium statistical physics, 
we expect scale invariance only in presence of certain 
symmetries or at critical points We are therefore lead 
to seek for the origin of the scale invariance in nature in 
the rich domain of nonequilibrium systems ■ One 

might hope, in fact, that there are classes of nonequilib- 
rium systems that generate scale invariance for a wide 
(and arbitrary) range of physical parameter, providing 
an explanation for the commonly observed scaling laws. 

Pursuing this aim, Bak, Tang and Wiesenfeld proposed 
the concept of Self-Organized Criticality (SOC) Q as a 
unifying framework to describe a vast class of dynam- 
ically driven systems which evolve spontaneously in a 
stationary state with a broad power law distribution of 
discrete energy dissipating events. To illustrate the ba- 
sic ideas of SOC, they introduced a cellular automaton 
model of sandpiles. In this model, criticality seems to 
emerge automatically if the system is driven at an in- 
finitesimal rate Because of the enormous concep- 
tual potentiality, SOC ideas have reverberated rapidly 
throughout the sciences, from geophysics to economics 
and biology, as a prototype mechanism to understand 
the manifestation of scale invariance and complexity in 
natural phenomena. 

The major source of difficulties in the study of sand- 
piles models is the absence of a general criterion, like 
the use of the Gibbs distribution in equilibrium systems, 
to assign an ensemble statistical measure to a particular 
configuration of the system. This problem is common to 
many nonequilibrium systems whose theoretical under- 
standing lies far behind the equilibrium theory. In partic- 
ular, many relations among sandpiles automata and sys- 
tems with nonequilibrium absorbing critical point have 



been recently enlightened jlO| . 

In the past years, we developed a renormalization 
group (RG) strategy for sandpile models 0] which has 
also been applied |Q to forest-fire models |^||. This 
approach deals with the critical properties of the sys- 
tem by introducing in the renormalization equations a 
dynamical steady state condition which provides non- 
equilibrium stationary statistical weights to be used in 
the calculation. This scheme, named the Dynamically 
Driven Renormalization Group (DDRG) 0], has been 
successively generalized as a renormalization framework 
for systems with a non-equilibrium critical steady-state. 
Recently, the DDRG approach has been improved includ- 
ing higher order proliferations through a general scheme 
[ ]l4| |. The method has also been applied to one dimen- 
sional sandpiles fl5|| and other non-equilibrium systems 

Here we discuss the application of the DDRG to sand- 
pile models, deriving systematically the previous RG 
schemes and presenting new results. We will in- 

troduce the general strategy of the DDRG for the critical 
height sandpile models, and its practical implementation 
for increasingly complex proliferation schemes. In order 
to treat such a high level of calculation complexity, we 
introduce a generating function for the basic recursion 
relations. The scheme is then extended by exploiting the 
analogy with a particular chain chemical reaction. Fi- 
nally its application to the calculation of the boundary 
critical behavior is shown. 

The paper is organized as follows: in Sec. 2 we intro- 
duce the class of sandpile automata and its mapping into 
a general nonequilibrium cellular automaton (CA). Sec. 3 
presents the Dynamically Driven Renormalization Group 
general scheme. Sec. 4 shows the explicit application of 
the DDRG to the sandpile in its simple scheme. In Sec. 
5 we present the actual calculations of the renormaliza- 
tion equations and their generating function and results 
obtained. Sec. 6 describes the extended chemical reaction 
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scheme and its results. Sec. 7 is devoted to the renormal- 
ization analysis of the boundary critical behavior. Sec. 8 
presents the summary and conclusions. 

II. THE SANDPILE MODEL 

The prototype example for SOC is provided by sand- 
piles: sand is added grain by grain until unstable sand 
(too large local slope of the pile) slides off. In this way 
the pile reaches a steady-state, in which additional sand 
grains fall off the pile by avalanche events. The steady- 
state is critical since avalanches of any size are observed. 
This class of models can be used to describe a generic 
avalanche phenomenon, interpreting the sand as energy, 
mechanical stress or heat memory. 

Sandpile models are cellular automata |t],|17| defined in 
a d— dimensional lattice. A discrete or continuous vari- 
able E(i), that we denote by energy, is associated with 
each lattice site i. At each time step an input energy 6E 
is added to a randomly chosen site. When the energy 
on a site reaches a threshold value E c , the site relaxes 
transferring energy to the neighboring sites 

E(i)^E({)-J2^(e) (1) 



N 

(s\T\s°)=l[r(8 i \s i ,{a i+e }) (3) 
»=i 

where N — L d is the number of sites, and e specifies the 
nearest neighbor (n.n.) vector. The dynamics is therefore 
expressed as a product of one-site transition probabilities 
depending upon the sites and its nearest-neighbors states 
at the previous time step. 

As we said, the common characteristic of SOC systems 
is the presence of a nonequilibrium critical steady-state, 
which we can analyze using the DDRG formalism. How- 
ever, it is worth remarking that SOC systems reach true 
criticality just in the limit of an infinite slow driving con- 
dition. This means that the perturbing time scale is much 
larger than the dynamical activity one. SOC systems re- 
lax far more rapidly than they are perturbed. In practice, 
this implies that no new grain of sand is dropped until 
the avalanches started by the previous grain has finished. 
In this way avalanches cannot overlap, and their dynam- 
ics is well defined with respect to the external field. A 
complete RG analysis should take into account also the 
driving field. However, since we are interested in the crit- 
ical point, we will study the system in the limit of slow 
driving. A more detailed discussion of the complete sand- 
pile automaton phase diagram is provided in ref. [fLO] . 



E(i + e) -> E(i + e) + AE(e) (2) 

where e represents the unit vectors on the lattice. A typi- 
cal choice for the parameters is , for example, E c = 4 and 
AE(e) = SE = 1, but other possibilities have also been 
considered. The relaxation of the first site can induce a 
series of relaxations generating an avalanche. Note that 
the energy is added to the system only when the configu- 
ration is stable (i.e. all the sites are below the threshold). 
The boundary conditions are usually chosen to be open 
so that energy can leave the system. In these conditions 
the system organizes itself into a stationary state charac- 
terized by avalanches of all length scales. In particular, 
the distribution for avalanches sizes s decays as a power 
law P(s) ~ s~ r , and the linear size of the avalanche 
scales with time r ~ t z . This model has been extensively 
studied in the past by means of numerical simulations 
]T8| pH and several exact results have been derived for 
abelian sandpiles models (ASM) p2|-p6[. 

Given the above definition of sandpiles we can rephrase 
them in the language of discrete nonequilibrium proba- 
bilistic CA. To each site i is associated a variable s$ that 
can assume q different values (s$ = 1, 2, 3, • • • , q). For in- 
stance, each state might correspond to an allowed energy 
level. The subscript i labels the lattice site. A complete 
set s = {s^ of lattice variables specifies a configuration 
of the system. We define (s | T(/i) | s°) as the transition 
rate from a configuration so to a configuration s in a time 
step t as a function of a set of parameters /j, = {fii}- SOC 
automata are usually defined by a transition probability 
given by the product 



III. THE DDRG 

The probability distribution of CA such as those shown 
in the previous section obeys the following master equa- 
tion (ME) 

P(s, h + t) = (s I T(p) | s Q )P( S °, t Q ). (4) 
W'} 

The explicit solution of the master equation is not in 
general available but we can extract the critical prop- 
erties of the model by a renormalization group analy- 
sis. We coarse grain the system by rescaling lengths 
and time according to the transformation x — > b x and 
t — * b~ z t. The renormalization transformation is con- 
structed through the operator 1Z(S, s) that introduces a 
set of coarse grained variables S = {Si} and rescales the 
lengths of the system Q. In general, TZ is a projec- 
tion operator with the properties TZ(S, s) > for any 
{Si},{si}, and X]{S}^(^ S ) = 1- These properties pre- 
serve the normalization condition of the renormalized dis- 
tribution. The explicit form of the operator TZ is defined 
case by case in various applications of the method. Usu- 
ally, it corresponds to a block transformation in which 
lattice sites are grouped together in a super-site that de- 
fines the renormalized variables Si by means of a majority 
or a spanning rule. 

We subdivide the time step in intervals of the unitary 
time scale (to = 0) obtaining the coarse graining of the 
system as follows: 
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P'(S,t') = 

E w R(^)E {s . } (s I T h \n) | s °)p( s °,o) (5) 

where we have included the application of the opera- 
tor K and t' = b z t. The meaning of (s | T b * (fi) \ s°) 
has to be defined explicitly: the simplest possibility is 
b z = N where N is an integer number, and T N denotes 
the application of the dynamical operator N times. In 



general, since we are dealing with a discrete time evo- 
lution, we have to consider T b as a convolution over 
different paths, chosen by an appropriate condition. The 
detailed definition of the effective operator T b ~ for the 
sandpile is reported in the next section. By multiply- 
ing and dividing each term of eq. ((H) by P'(S°,0) — 
J2{s°} 7Z(S°, s°)P(s°, 0) and using the properties of the 
operator TZ, after some algebra we get 



p'(s,t') = £(- 

{s°} 



E {s o } E M k{S°, 8°)n{s, \ t* z m | s °)p( s °, o) 



E {s c } K(5VW,0) 



)P'(S°,0) 



(6) 



which finally identifies the renormalized dynamical oper- 
ator (S | T' | S°). In other words, the new dynamical 
operator T' is the sum over all the dynamical paths of 
b z steps that from a starting configuration {s^} lead to 
a configuration {si} which renormalizes respectively in 
{Sf} and {Si}. The sum is weighted by the normalized 
statistical distribution of each configuration. The scheme 
discussed so far is a general formulation valid for each 
system which exhibits a stationary state, and its appli- 
cation presupposes the knowledge of the explicit form of 
the steady-state distribution W(s) = P(s,t — » oo). For 
instance, in equilibrium phenomena W(s) is given by the 
Gibbs distribution. In this case it is possible to apply sev- 
eral methods such as cumulant expansions and exact or 
approximate decimation to obtain the form of the recur- 
sion relations. For non-equilibrium dynamical systems, 
in general we do not know the form of the steady-state 
distribution. We will therefore develop an approximate 
method to evaluate the stationary distribution to be used 
in the calculation of the renormalized master equation. 

The steady-state distribution can in general be split 
into two parts 

(c) 



(7) 



W(s) = W®(s) + WW{s) 

where W^\s) and w( c \s) are, respectively, the inco- 
herent and coherent part of the distribution. The inco- 
herent part of the distribution does not include correla- 
tions among variables and expresses a mean field approx- 
imation for the system. The coherent part VF( c )(s) can 



be subdivided in parts describing different kinds of cor- 
relations: nearest-neighbors, next-nearest-neighbors, etc. 
The incoherent part is a factorized distribution that, for 
systems characterized by a q-state variables (see sect. 2), 
has the form 



where (p K ) is the average density of sites in the K-state. 
In this way, we have approximated the probability of each 
configuration {si} as a product measure of the mean field 
probability to have a state Si in each corresponding site. 
The incoherent part contribution to the renormalization 
equation can be obtained by stationarity conditions for 
the system ^({(/Ok)}) = to evaluate the densities (p K ). 
These conditions are derived from dynamical mean field 
equations which describe the driving of the system to the 
nonequilibrium steady state by means of balance con- 
straints. The operator depends upon the same dy- 
namical parameters of the operator T, and by solving 
the stationary condition equation, the average densities 
of the K-states for the coarse grained system are obtained 
as a function of /i at the corresponding iteration of the 
RG equations. By inserting the approximate distribution 
in Eq. (^), we thus get the following set of renormaliza- 
tion equations 



(S | T'(/i) | S°) 



E {s o } ^(5o, s o)n,(p s o) 



(9) 



SM( Pk )}) = 



(10) 



where the second equation denotes the dynamical steady 
state condition that allows evaluation of the approximate 
stationary distribution at each coarse graining scale. We 



call Eq. ( pp| ) the driving condition, since it drives the 
RG equations acting as a feedback on the scale trans- 
formation. Eqs (||) and (|l^) are the basic renormaliza- 
tion equations from which the desired recursion relations 
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are derived. Imposing that the renormalized operator 
T' has the same functional form of the operator T, i.e. 
T'(p) = T(p'), we obtain the rescaled parameter set 
p! = f(fJ>). This implies that the renormalized single 
time distribution P'(S,t') has the same functional form 
of the original distribution P(s,t). The critical behav- 
ior of the model is obtained by studying the fixed points 
A 4 * = /(A 4 *)- Since we are dealing with discrete evolu- 
tion operators T, we define the time scaling factor b z 
as the average number of steps and apply the operator 
T in order to obtain that T'(p) = T(p') for the coarse 
grained system. In this way we obtain a time recursion 
relation t' = g(p)t, or equivalently b z = g(p), from which 
it is possible to calculate the dynamical critical exponent 
z = log g(p*) /log b. In this form of the DDRG, we take 
into account only the uncorrelated part of the steady- 
state probability distribution. The results obtained are 
not trivial because correlations in the systems are con- 
sidered in the dynamical renormalization of the operator 
T, that given a starting configuration traces all the possi- 
ble paths leading to the renormalized final configuration. 
Moreover, geometrical correlations are treated by the op- 
erator 1Z that maps the system by means of spanning con- 
ditions or majority rules. The renormalized uncorrelated 
part of the stationary distribution is evaluated from the 
stationary condition with renormalized parameters, thus 
providing an effective treatment of correlations. One can 
then improve the results by including higher order con- 
tributions to the unknown stationary distribution W(s) 
using cluster variation methods 



IV. RENORMALIZATION SCHEME FOR 
SANDPILE MODELS 



Here, we show in detail the DDRG scheme for sand- 
pile models. For the sake of clarity, we start by consid- 
ering the minimum proliferation scheme. A more refined 
scheme is discussed in the following sections. 

To simplify the description of sandpile models as much 
as possible, we can reduce the number of states of each 
site in the following way. At any scale, we divide the 
sites in critical (sj = 1) and stable (s$ — 0). Stable 
sites do not relax when energy is added to them. On the 
other hand, critical sites relax when they receive an en- 
ergy grain 5Ei n . In this formalism we define < p > as the 
density of critical sites. For convenience, we also define 
unstable sites (s, = 2) , as those that are relaxing, even 
though they are not present in the static configurations 
of the system. These definitions can be extended to a 
generic scale b. For instance, a cell at scale b is consid- 
ered critical if the addition of energy 5Ei n (b) induces a 
relaxation of the size of the cell (i.e. the avalanche spans 
the cell). 

In a relaxation event at the minimal scale energy is 
equally distributed in the four directions. This is no 
longer the case at a coarse grained level where differ- 



ent possibilities arise: the energy in principle can be dis- 
tributed to one, two, three or four neighbors. It is also 
worth remarking that in a certain case unstable sites at 
the coarse grained scale do not dissipate energy to nearest 
neighbors, representing just intra-site energy rearrange- 
ments. These processes define the probability that relax- 
ation events take place on the renormalized scale without 
energy transfer. All these events occur with probabilities 

P = (P0,Pl,P2,P3,P4) 

In terms of the matrix element (0|T|2) the vector P rep- 
resents the probabilities 



Pn 



\T\2)r 



(11) 



where (0|T|2)„ is the probability that a relaxing site be- 
comes stable and transfers energy to n neighbors. In 
this way, we have obtained the set of parameters that 
describes the dynamics. Of course, the choice of the pa- 
rameter space is not uniquely determined; one encounters 
proliferation problems typical of real space RG methods. 
For instance, higher order proliferations are due to mul- 
tiple relaxations of the same site and sites becoming crit- 
ical during the dynamical process (i.e.: (1|T|2)). In the 
following, the practical implementation of the method 
considers just the minimal proliferation we have reported 
above. In the next sections, a more refined scheme will 
be treated. 

First of all, let us show how the driving condition is ob- 
tained by imposing the stationarity of the process. The 
average energy of a site evolves according to the following 
equation written in the continuum notation: 



dE(t) 
dt 



8E„ 



5E n 



(12) 



where SE in is the average energy entering into the site 
either because of relaxation in a neighboring site or be- 
cause of the external perturbations, and SE out is the av- 
erage energy dissipated by the site. The stationary state 
is characterized by the balance between the energy that 
goes in and the energy that goes out of the system. Wc 
assume that energy is transferred in "quanta" 5E = 5Ei n 
in each direction and we on average obtain 



which implies 



5E={ P )5EJ2 



(P) 



np n 



T,n n Pn 



(13) 



(14) 



This relation gives the average density of critical sites in 
the steady-state, allowing us to evaluate the approximate 
stationary distribution at each scale. 

The renormalized matrix element is then obtained by 
considering all the renormalized processes that span the 
cell and transfer energy outside 
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Pn = {Si = 0\ T' | Sf = 2) n . (15) 

We proceed in denning explicitly a renormalization pro- 
cedure for the dynamics by considering a finite trunca- 
tion on four-sites cells. This corresponds to a cell-to-site 
transformation on a square lattice, in which each cell at 
the coarser scale is formed by four sub-cells at the finer 
scale: the length scaling factor is b — 2. In this case, the 
operator 1Z can be written in the following way: 

TZ(S,s) = l[n(Sj,{s i }j) (16) 
J 

where each term is acting on a specific cell J and {si}j 
denotes the configurations of sites belonging to that cell. 
A cell is renormalized as a relaxing one if it contains a 
relaxing sub-cell which transfers energy to a critical sub- 
cell. In this way, we ensure that the occurring relaxation 
process is extending over the size of the renormalized 
length scale independently of the successive avalanche 
evolution. A critical cell is therefore defined by a cell 
which can be spanned by a path of relaxation events. 
The scheme considers only connected paths that span 
the cell from left to right or top to bottom. This span- 
ning rule implies that only paths extending over the size 
of the resulting length scale contribute to the renormal- 
ized dynamics, and it ensures the connectivity properties 
of the avalanche in the renormalization procedure. 

Every cell at the coarser scale can be characterized by 
an index a that indicates the configuration of sub-cells, 
and we have that 53 { s } ~* S Q - The approximated sta- 
tionary distribution (Eq. (||)) for each of these configura- 
tions is given by 

4 

W a ((p)) = n a H(p Si ) (17) 

where n a is a factor due to the multiplicity of each con- 
figuration. 

By using this scheme and replacing sums over config- 
urations with sums over the index a, the recursion rela- 
tions can then be rewritten in the simpler form 

Pn = ^£^«P»£<«' I T& >«') I «>» ( 18 ) 
a a' 

where | a), \ a') denotes the four site configurations 
which renormalize in | Sf — 2) and | Si =0), respec- 
tively. In the above expression the denominator of eq. (j^) 
is adsorbed in the normalization factor J\f. 

The effective operator T b contains all the dynamical 
processes that contribute to the definition of a mean- 
ingful renormalized dynamics. We define the following 
transformation 

MT^I^H^P^IT^)!* } (19) 

N 

where T>n is the renormalization operator for the dynam- 
ical evolution of the system: it is a projection operator 



that samples only the paths of N time steps which have 
to be considered in the definition of the effective operator 

The operator T>n is chosen on the basis of physical 
considerations: spanning conditions etc. In addition, 
T> n should satisfy some general properties in order to 
preserve the symmetry or the internal space of the dy- 
namical variables. For instance, we have to ensure the 
normalization of the effective dynamical operator by the 
property 

££2> J vHT"(,x)| S °> = l (20) 

Moreover, T>n must be consistent with the definition of 
the renormalization operator 1Z: it should describe dy- 
namical processes among renormalized variables of the 
same type of those given by the operator 1Z. Finally, T>n 
has to preserve the form of the dynamical operator T at 
each scale. This condition imposes that the time scaling 
is consistent with the length scaling used in JZ. In this 
way, it is possible to map the renormalized system in the 
old one with renormalized variables. The operator T>jq 
is therefore defined explicitly as an operator acting on 
the paths internal to four site cells. It selects for each 
N just relaxation paths which consist of N connected 
non-contemporary relaxation events that leave the cell 
without unstable sites. In a mathematical form it reads 
as 

N-l 4 

v N = n (i-wiie^- e ( 2i ) 

i£{a'} J=0 m=l ie{aj} 

where aj's are the intermediate cell configurations dur- 
ing the dynamical evolution and J2ie{ aj } denotes the 
sum over all the sites in the cells. In the above expres- 
sion, each delta function acts on a different intermediate 
cell eliminating those paths which do not have activity at 
each dynamical step. Furthermore, the operator ensures 
that in the cell a' (Nth step) no activity is present; i.e 
the process has stopped. Finally, we have to write the 
equation that gives the time scaling factor from the total 
average over contributing processes to the renormalized 
matrix element (0 | T | 2) 

0(p») = ^£^«p»££^2W i i «> 

a a' N 

(22) 

where we used the DDRG scheme to explicitly get the 
stationary weights, and J\f is an opportune normaliza- 
tion factor. 

The above relations will provide the consistent rescal- 
ing of time by imposing that b z = g(Pn) from which it is 
possible to calculate the dynamical critical exponent. 
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V. MINIMAL PROLIFERATION 
CALCULATIONS AND RESULTS 

The explicit evaluation of the recursion relations de- 
pends on the choice of the spanning condition. In the fol- 
lowing, the scheme used considers only connected paths 
that span the cell from left to right or from top to bottom. 
This spanning rule implies that only paths extending over 
the size of coarse grained length scale contribute to the 
renormalized dynamics, and it ensures the connectivity 
of the avalanche in the renormalization procedure. 

An example of such a path is shown in Fig. (|l|) . In this 



case a = 2 and the path shown refers to the probabil- 
ity that the unstable sub-cell relaxes towards the other 
critical sub-cell [Fig. 1(b)]. This occurs with probability 
(l/4)pi. At this point we consider the probability that 
the next relaxation event at the fine scale involves two 
neighboring sites, one inside and the other outside the 
original cell of size b — 2 [Fig. 1(c)]. This occurs with prob- 
ability (2/3)p 2 . This series of relaxation processes con- 
tributes to the term V 2 {a' | T 2 (p n >) \ a = 2)i that char- 
acterizes the relaxation processes at the coarse grained 
scale. By summing over all the paths that lead to a p[ 
process one obtains for a = 2 



J2a'J2N^N(a'\T N (p nl )\a^2) 1 = 

(jPl + \P2){\pi + fP2lP3) + (fP2 + \pi)(\pi + \P2)+ (23) 
{\P2 + lP3)(|j>a + jPs)(|Pl + \P2JP3,)- 



In a similar way one can also write the expression for the 
complete recursion relations. Proliferation effects due to 
multiple relaxations of the same site and sites becoming 
critical during the dynamical process are not considered 
in this scheme. However, the complete polynomials for 
p' n involve more than two hundred terms, and we devel- 
oped a generating function that allows their systematic 
calculation. 

The generating function allows us to find the form of 
the renormalized operator T', without writing down the 
explicit form of all relaxation paths in the coarse grained 
cell. The basic idea of this method is to renormalizc 
the function which describes all relaxation paths at once, 
rather then the probabilities of separate relaxation paths, 
by using the branching structure of avalanches, 

To describe in detail the branching process underly- 
ing the large scale behavior of the sandpile model, we 
consider the generating function 

a(N, E,S,W) =p + ?j(N + E + S + W) 

+ ^-(NE + NS + NW + ES + EW + SW) (24) 
6 

+ ?j{NES + NEW + NSW + ESW) +p 4 NESW, 

where symbols N, E, S and W correspond to the north, 
east, south and west directions on the square lattice, re- 
spectively. The coefficient in each term of this polynomial 
gives the probability for the process to go in the corre- 
sponding directions. The generating function takes into 
account all possible relaxation processes in the cell. It is 
easy to check directly that this function has the following 
properties: 

1. If the argument corresponding to any direction is 
replaced by zero, the function counts the relaxation 
processes that do not send energy to this direction 
(Fig.|(b)). 



2. If the argument corresponding to any direction is 
replaced by unit, the function counts the relaxation 
processes whether or not the energy is transfered in 
that direction (Fig.S(c)). 

3. The generating function is normalized so that 
<t(1,1,1,1) = 1. 

If there is a critical cell near the relaxing one, the 
outgoing energy can initiate the relaxation of the cell. 
It is easy to see that we can replace the argument 
corresponding to this direction by another generating 
function corresponding to the relaxation of the second 
cell. Finally, we obtain the generating function of this 
two-step relaxation process. For example, the func- 
tion <r(JVi, (r{N 2 , E 2 , S 2 , 1), Si, Wi) describes the pro- 
cesses where the cell 1 relaxes first. Then, if the energy 
goes to the direction N\ it initiates the relaxation of the 
cell 2 (fig.|(d)). The symbols N l ,E l ,S l ,W l denote the 
directions outgoing from the cell i. 

Using these properties we can write down the generat- 
ing function £„, counting the relaxation processes in the 
block that consists of four cells for different a. To this 
end, we must take into account only the processes which 
match the spanning condition. Therefore, it is necessary 
to eliminate all processes, in which only one cell relaxes. 
As the coefficients of the polynomial E a have the mean- 
ing of probabilities, they should finally be normalized by 
the condition 

S Q (1,...,1) = 1. (25) 

The generating function corresponding to the relax- 
ation processes inside the block with a = 2 is 

S 2 = 

{a{a{N 2l 1, 1, W 2 ), 1, Si, Wi) - <r(0, 1, Si, Wi)} + 
{a{N 2l l,a(l, 1, Si, Wi), W 2 ) - a(N 2 , 1, 0, W 2 )} + (26) 



G 



c.p.)/Z 2 , 

where Z 2 is a normalization factor chosen so that 
£2(1, •••! 1) = 1- To write this function we start from 
the left down cell and define the arguments of the a- 
function corresponding to the toppling of this cell. For 
the process definitely spans the block, the left up cell 
should topple and we write another u-function instead of 
the symbol N±. By going eastward the process will ter- 
minate inside the block and this branch of the toppling 
process cannot affect the neighboring blocks. Hence, we 
should write the number 1 instead of the symbol E\ . The 
other symbols Si and W\ correspond to the branches of 
the relaxation process that goes immediately out of the 
initial block of cells. To consider only the processes that 
span the block, we must subtract the cr-function describ- 
ing processes that do not send the energy from the first 
critical cell to the second one. Then, we add analogous 
cr-functions for the processes starting from the relaxation 
of the left upper cell. The term c.p. denotes all possible 
cyclic permutations of the critical cells inside the block. 
Analogously, we can write cr-functions of all other types 
of blocks. 

To obtain complete generating function £ for the block 
of four cells, we should sum up all £ Q -functions with the 
weights of blocks and normalize the result 

£ = |]T^ Q £ Q . (27) 

OL 

Now, to transform the ^-function from the block of four 
cells at the scale b k to larger cell at the next scale b k+1 , 
we replace the directions N\,N2,... outgoing from the 
initial block by the new arguments corresponding to the 
directions N, ... outgoing from the new renormalized cell. 
In other words, two bonds that connect the neighboring 
blocks are coupled to the only bond on the lattice at the 
next scale, as is shown in fig]|. Eventually, we obtain the 
following generating function: 

E(JV, E, S, W) = 2 2 ^ 3 3+VV4 4 (28) 

Zi 

where 

£ 2 = 

(a(a(N, 1, 1, W), 1, S, W) - (7(0, 1, S, W)) + 
(a(N, 1, (7(1, 1, S, W), W) - a(N, 1, 0, W)) + c.p.)/Z 2 , 

(29) 

X 3 (N,E,S,W) = 

(a(a(N, a{N, E, 1, 1), 1, W), 1, S, W) - a(0, 1, S, W) 
+ a{a{N, 1, 1, W), <r(l, E, S, 1), S, W) - ct(0, 0, S, W) 
+ ct(1, a(a(N, E, 1, 1), E, S, 1), S, W) - a(l, 0, S, W) 
+ c.p.)/Z 3 

E 4 (N,E,S,W) = 

(a(a(N, a(N, E, a(l, E, S, 1), 1), 1, W),0, S, W) + 



cr(0, a(a(N, E, 1, a(N, 1, 1, W)), E, S, 1), S, W) - 

2a(0, 0, S, W) + (d(l, 1, S, W) - a{l, 0, S, W) - 

<t(0,1,S, W) + (7(0,0,5, W)) 

(<j(N, 0, 1, W)a(0, E, S, 1) + a(N, E, 1, 1) 

(<j(N, 0, 1, W)(a(l, E, S, 1) - (7(0, E, S, 1)) + 

(7(0, E, S, l)(a(N, 1, 1, W) - a(N, 0, 1, W)) + 

(a(l,E, S,l)-a(0,E, S,l)) 

(a(N, 1, 1, W) - a(N, 0, 1, W)))) + c.p.)/Z A 

Here, Zi and Z are the normalization factors and c.p. 
denote the expressions obtained from the previous poly- 
nomial by all possible cyclic permutations of its argu- 
ments. This generating function is the polynomial that 
contains only the first and second powers of its argu- 
ments. The last terms correspond to the processes when 
two energy portions are transferred from the initial block 
to the neighboring block by the two paths. However, ac- 
cording to the RG strategy, these processes should be 
considered as the transfer of the coarse grained energy 
portion at the larger scale. Therefore, all second powers 
of the arguments should be replaced by the first ones. 
The result obtained is the generating function describing 
the relaxation of renormalized cells. It depends on the 
same products of its arguments as the generating func- 
tion for initial cells, but the coefficients of this polynomial 
are different and equal to the probabilities of relaxation 
processes that send energy to the given direction at the 
new scale. Taking these coefficients, we obtain the sought 
recursion relations which link the parameters of the cell 
at the scale b k with the same parameters at the scale b k+1 

P{k + l) = f{P{k),p(k)) (30) 

The above set of RG equations supplemented with the 
driving condition 

(P k+1 ) = ^-TTT ( 31 ) 

define completely the DDRG recursion relations for sand- 
pile models. Given this scheme the flow diagram and the 
relative fixed point in the parameters space (p, P) can be 
studied. We consider here the calculation scheme imple- 
mented with P = (j>o,pi,p2,P3,P4 : ). Despite the enlarge- 
ment of the phase space by including the proliferation 
characterized by the probability po the flow in the phase 
space is very similar to those obtained in ref. JHJ, where 
this parameter was not considered. A single attractive 
fixed point is obtained and the numerical value of this 
fixed point is very close to that obtained in the approach 
of Ref. |ll|. The complete attractiveness of the fixed 
point corresponds to the lack of relevant scaling field, i.e. 
control parameter. This must be the case, because we 
implement our RG scheme in the infinite time scale sep- 
aration limit. This implies the relevant scaling field are 
constrained to their critical values. In other words we are 
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restricting the study of the system on its critical surface. 
In Tab. I we report the results obtained within this cal- 
culation scheme. The single fixed point is the signature 
of a single universality for all the non-directed sandpilc 
models. This is questioned from numerical simulations 
which show evidences this could not be the case |3^| . It 
appears, in fact, that BTW model and Manna model 
belong to different universality classes. This distinction 
does not appear in the present proliferation scheme. The 
present approach, however, can not provide a definitive 
settling of this issue. Mainly, it depends on the fact that 
we still are neglecting some proliferations such as the pos- 
sibility of multiple topplings, that could be relevant for 
the identification of different universality classes. 

The avalanche exponent t can be obtained directly 
from the fixed point parameters. By using the discrete 
length scale = 2 k and the avalanche distribution in 
the form P(r)dr w A x ^ 2T ^dr we can define the probabil- 
ity that the relaxation process spans the cell of size 
and dies at the neighboring cells not extending over the 
scale M fc+1 ) 



I P(r)dr/ / P{r)dr = 1 - 2 2(1 ~ t) (32) 

(,(*-!) Jb( k - 1 ) 



K = 



Asymptotically (k — > oo) we can express K in terms of 
fixed point parameters p* and p* in the following way: 



K=p*+p* 1 (l-p*)+P* 2 (l-p*f 
+ P* 3 (l-P*f+Pl(l-P*) i . 



(33) 



This equation gives the total probability that a relax- 
ation process occurs without triggering other sites, and 
therefore it does not extend on length scales larger than 
that of a single cell. Using these two expressions, Eqs. 
(p4p3|), the exponent r is given by the formula 



In 1 -K) 

t = 1 i '- = 1.262 

2 In 2 



(34) 



in excellent agreement with the proposed value t = 5/4 
p8[ and large scale numerical simulations p8| , ^9|j2l| . The 
obtained value is also in good agreement with the value 
obtained in the calculation of ref |ll[ , showing the robust- 
ness of the method with respect to different proliferation 
scheme. In order to overcome some of the approximations 
considered so far, we will present in the next section an 
improved scheme which takes into account a wider set of 
dynamical parameters. This scheme allows us to study 
also the critical behavior at the boundary of the system. 



VI. EXTENDED KINETIC EQUATION SCHEME 

In this section we treat more explicitly the dynamics 
of the original sandpile model. Considering the evolution 
of stable cells, we can take into account some of the pro- 
cesses which were neglected in the scheme discussed in 
the previous section. 



To keep the connection with the original formulation of 
the sandpile model, we will characterize the static prop- 
erties of a cell by four quantities 

N(k) = (n A , n B , n c , n D ), n A + n B + n c + no = 1, 

(35) 

which are nothing but the probabilities for a cell to be- 
have like a site on the initial lattice with a height 1, 2, 
3 or 4, respectively, in the coarse grained dynamics, i.e. 
the addition of a "coarse grained particle" to the cell 
transforms it to the next one in the alphabet. For exam- 
ple, the cell B characterized by the vector (0,1,0,0) will 
be transformed to the cell C with the vector (0,0,1,0). 
The last variable «d is the probability for the cell to be- 
have like a critical one in a sense that the addition of a 
"coarse grained particle" to the cell induces relaxations 
into some neighboring cells or, in other words, subre- 
laxation processes on a minimal scale span the cell and 
transfer energy to some of its neighbors. 

As we stressed in the previous sections, independently 
of the dynamics of the model at the minimal scale, each 
critical cell is characterized by the vector 

P( k ) = (Pl,P2,P3,Pi), Pi +P2+P3+P4 = 1, (36) 

that gives the probabilities for the energy to go to 1,2,3 
or 4 neighboring cells after the relaxation of the critical 
cell. Here, because we have already enlarged the phase 
space by introducing the densities n a , we do not include 
Po in the calculation scheme. 

In this framework the coarse grained dynamics of 
the sandpile model can be represented as the following 
branching process on the sublattice Lb'. 



A 
B 



C + p 



(37) 



D 
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Pi ■ 


D- 


h£ 


P2 ■ 


C- 


-2p 


P3 ■ 


BH 


- 2>(p 


Pi ■ 


A- 


h4(^ 



Here, (p and dp denote the "coarse grained particles" ob- 
tained by the cell and the particles transferred to the 
neighboring cells, respectively. 

These processes can be formally reinterpreted as an 
irreversible chemical reaction which takes place at each 
cell of the sublattice Lb- Now the coarse grained variables 
nAjiB.^c^D and n v denote the concentrations of the 
respective species A, B, C, D, and ip. Following to stan- 
dard prescriptions of the chemical physics we can write 
down kinetic equations corresponding to this scheme of 
chemical reactions 
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"A 


= n v 


(jpi n D 


- "A 


), 


(38) 


flB 


= n v 


(P3 nr> - 


f n A 


- n B j, 


(39) 


h c 


= n v 


(P2 n u - 


f n B 


- n c ), 


(40) 


hu 


= n v 


(pi n D - 


f n c 


~ «d), 


(41) 


n v 


= n v 


(p n B - 


1) + 


p v\7 2 (n v n D ) +r)(r,t) 


(42) 



where p = p\ + 2j>2 + 3j>3 + 4j»4 is equal to the average 
number of particles leaving the cell on toppling and r 
is the position vector of the cell in the 2D space. The 
noise term t](r, t), being non- negative, mimics the ran- 
dom addition of particles to the system. The diffusion 
term V 2 (n ¥ ,riD) describes the transfer of particles into 
the neighboring cells, and the diffusion coefficient v for 
the discrete Laplacian on the square lattice is equal to 
1/4. 

The only mobile specie in this scheme of reactions is 
ip and it is the field n v which describes the dynamics of 
avalanches. When it is equal to zero, all toppling pro- 
cesses die. Then, due to the noise term r](r,t), particles 
are added randomly into the system initiating a branch- 
ing process directed to the open boundary of the system. 
This process mutates species in the cells it has visited 
and topples the critical ones. Finally, the system will 
reach the steady state where the probability that the ac- 
tivity will die is on average balanced by the probability 
that the activity will branch. Thus, the chain reaction 
maintains this stationary state and all further avalanches 
cannot change the concentrations of species A, B, C, and 
D. Therefore, the steady state is characterized by the 
conditions that 

n A = n B = n c = n D = (43) 

and Eqs. ( [38||4T| ) lead to the following relationships be- 
tween concentrations of species N(k) at the stationary 
state and branching probabilities P(k) 



n *A = Pi/P, 

n% = {P3+Pi)/P, 

n C = {P2+P3 +Pa)/P, 



It 



D 



(Pi +P2 +P3+Pi)/P = 1/P 



(44) 

(45) 
(46) 
(47) 



The relation ( f47j ) between the probability and branch- 
ing probabilities P{k) can also be derived from the as- 
sumption that at the stationary state the flow of particles 
in a cell was on average balanced by the flow of particles 
out of the cell. 

Thus, we have found the driving conditions for the 
sandpile models. Using them we can link the statistic 
weights of any static configurations of a cell with the 
dynamic parameters. Now, we can realize the renormal- 
ization procedure described previously. To this end, we 
must consider all types of blocks of four cells, whose re- 
laxation matches the spanning condition. Such blocks 
and some of their relaxation schemes are shown in fig^. 
While the previous scheme deals only with the cells be- 
ing critical before the relaxation of the block, this one 



allows us to consider the cells becoming critical during 
the relaxation(fig.||(c)). 

To obtain the recursion relations by the method pre- 
sented in the previous section, it is necessary to calculate 
the statistic weights of all configurations considered. The 
statistic weight of the block is given by the product of 
probabilities rii for all cells in the block, multiplied by 
the number of different blocks with the same relaxation 
schemes. Thus, the following weights must be ascribed 
to the blocks shown in fig.f|: 



W a 
W h 
W c 
W d = r&. 



4no(n A 
4nl(n A 



n B ), 



(48) 
(49) 
(50) 
(51) 



Expressing them through the probabilities P, by using 
driving conditions ( p4[]47| ), we obtain the complete sys- 
tem of renormalization equations. 



P(k + 1) = /(P(fc)) 



(52) 



Given this set of RG transformations we can study how 
the system evolves under successive doubling of length 
scale. The final result is independent upon the initial 
conditions ni and pi at the minimal scale. Also in this 
case the fixed point is attractive in the whole phase space 
and the system evolves spontaneously toward the fixed 
point values p* and n* shown in Tables || and III . These 
results can be compared with the exact ones obtained for 
the sandpile j24|. The exact height probabilities for the 
sandpile are reported in Tab. Ill and compared with our 
RG results. In order to calculate the avalanche exponent 
we can use the eq.([32|) by expressing K in terms of the 
fixed point parameters in the following way: 



K=pl{l-nl)+p* 2 {\- 
+ pl(l-n* D ) 3 +pl(l 

More generally the probability (|53 
via CT-function 



-n* B )\ (53) 
) should be represented 



K = o*(l- n* D (N), 1 - n* D {E), 1 - n* D {S), 1 - n* D (W)) 

(54) 

where no(N), Ud{E), ri£,(S) and Ud(W) are the concen- 
trations of critical cells at the nearest neighbors. By us- 
ing the fixed point values we finally obtain r = 1.248, 
that again is in very good agreement wit h the pro posed 
value 5/4 



19 21fl. Also 



and numerical simulations 
in this case the scheme results in a single universality 
class for the sandpile models. In the extended scheme 
we include part of the proliferations by allowing different 
heights, and taking into account that some of the sites 
becomes critical during the relaxation event. We do not, 
however, consider multiple relaxation of the same sites 
during the spanning time, nor we allows a renormaliza- 
tion of the energy transfer (8E). This parameters could 
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be important in the case of the Manna model as pointed 
out in ref . |32j . Work is in progress to extend the present 
DDRG scheme in order to include also these further pro- 
liferations. 



VII. THE BOUNDARY CRITICAL PROPERTIES 

Since the critical properties of the sandpile model are 
quite similar to those of second order phase transitions, 
we proceed here along the same lines followed in the 
study of equilibrium critical phenomena. In particular 
we determine also the surface critical exponents, which 
in general differ from the bulk ones. This is of special 
importance in the two-dimensional case where conformal 
field theory connects surface and bulk properties of the 
model 131. 



A. Open boundary 

The fact that the boundary is open means that after 
the relaxation of the boundary site the energy can leave 
the system. We consider the critical energy of sites at the 
open boundary to be E c = 4. It is more convenient to 
consider the boundary lying along the diagonal of the lat- 
tice and construct the renormalized cell in the following 
way. We consider the block of cells 2x2 that contains one 
bulk, two boundary and one external cells as shown in fig. 
||. For the critical properties of the model at large scales 
does not depend on the local structure of the lattice, the 
results obtained should not depend on the specific choice 
of the boundary. To describe the boundary cells at an 
arbitrary scale, we introduce the vectors 



N + E 



(57) 



N°(k) = (n° A ,n° B ,n° c> n° D ), 

O i O i O i O i 

n A + n B +n c + n D = 1, 



Statistical weights of the static boundary configurations 
are given by the product of the probabilities of one bulk 
and two boundary cells. Now, the coefficients of the gen- 
erating function of the renormalized cell gives us the re- 
cursion relations 



P°(k + l)=f (P°(k),P(k)). 



(58) 



Together with the bulk recursion relations of Eq. (|52j), 
Eq. ([58]) represents the complete system of the renormal- 
ization equations for the case of the open boundary. This 
system also has only one fixed point. The obtained fixed 
point parameters are given in Tables |TJ and III. The 
comparison of fixed point height probabilities with the 
exact values obtained for ASM |25| shows rather good 
agreement. 

To calculate the boundary critical exponent we should 
use the generating function describing the relaxation of 
the boundary cell in the fixed-point of RG flow. It is 
given by eq. (32), where 



K = <(!,!, 1-7^,1-r^). 



(59) 



This indicates that the toppling of coarse grained bound- 
ary cell should stop in two neighboring internal cells. The 
result t = 1.486 is very close to t d = | calculated ex- 
actly in [^6| for open boundary of ASM. Such a good 
agreement is probably to be ascribed to the fact that in 
the boundary avalanche in ASM each site topples only 
once 1 27 . Therefore renormalization scheme that does 



(55) 



not take into account multiple topplings gives most real- 
istic results near open boundary. 



B. Closed boundary 



P°{k) = (pl^plpll 
P°i+P° 2 +P° 3 +Pl=l, (56) 

which have the same meaning as in the previous sec- 
tion. The kinetic equations and feedback relations coin- 
cide with the bulk ones. Thus, it is only necessary to 
find the correct form of recursion relations. To write the 
generating function for the boundary block, we can use 
again the generating functions describing the relaxations 
of bulk and boundary cells. Also we have to introduce 
a special generating function corresponding to the relax- 
ation of the unphysical external cell. Since there are not 
processes transferring energy from the external half-plane 
of the lattice to the internal one, we require that the ex- 
ternal cell immediately transfers the energy outside the 
lattice. Thus, we provide conservation of the flaw of en- 
ergy through the boundary in the scale transformation. 
For the block shown in fig.^j the generating function of 
relaxation of the external cell has the simple form 



Let us consider the sandpile model on the half-plane. 
The edge of this half-plain is the closed boundary di- 
rected along one of the lattice axes. Each boundary site 
has three neighboring sites. Two of them also belong 
to the boundary. The fact that the boundary is closed 
means that after the relaxation of a boundary site energy 
does not leave the lattice, being distributed among the 
neighboring sites. Since the energy can leave the bound- 
ary site only in three directions, it is quite natural to take 
critical energy of a boundary site E c = 3. 

In order to follow the RG strategy, we again perform 
the site-to-cell transformation, replacing the block of four 
cells by a single cell at the larger scale. The cells at an 
arbitrary scale can be considered as either boundary or 
bulk ones. While the former consist only of bulk cells 
at the smaller scale, the latter include both the bulk and 
the boundary cells of smaller size. The renormalization of 
the bulk cell is described by the system of RG equations 
obtained in the previous section. Let us introduce the 
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description for the dynamics of boundary cells at an ar- 
bitrary scale. The static states of a boundary cell can be 
represented by three symbols A,B,C, which correspond 
to energy values E = 1,2,3 of boundary sites on the ini- 
tial lattice. The addition of energy transforms the cell 
from the state A into the state B and the cell from the 
state B into the state C. The cell in the state C is critical. 
The addition of energy initiates its relaxation when the 
cell turns into the states A or B or remains in the state 
C, sending the energy to three, two or one neighboring 
cells, respectively. The probabilities for the cell on the 
closed boundary to be in one of the three states is given 
by the vector 

N c (k) = (n c A ,n B ,n c c ), n c A + n c B + n c c = 1. (60) 

In the same way, the probabilities for the energy to be 
transferred in one, two or three directions is given by the 
vector 



P c (k) = (p c 1 ,p c 2 ,p c 3 ), pl+p c 2 +p c 3 = l. 



(61) 



The relaxation process at the boundary cell can be rep- 
resented as follows: 



<r c (N, E, W) = ^(N + E + W) + 

^ (NE + NW + EW) + p%NEW. (70) 

o 

The general idea of the generating function for the block 
of four cells is the same as in the bulk case. The difference 
is that the generating functions for boundary blocks ori- 
ented differently with respect to the boundary should be 
calculated separately and cannot be obtained by simple 
cyclic permutations of the arguments. Finally, applying 
the renormalization procedure with the use of the feed- 
back relations (p6[pq), we obtain the recursion relations 



P c (k + 1) = f c (P c (k),P(k)), 



(71) 



where P(k) matches the bulk recursion relat ion s of 
Eq. ( p2| ) . The obtained height probabilities (Table III ) are 



in good agreement with those calculated exactly in the 
case of ASM with closed boundary. To calculate the crit- 
ical exponent r c describing the distribution of avalanches 
near a closed boundary, we use the Eq.(|3^), expressing K 
through the fixed point generating function of the bound- 
ary cell: 
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K = a* c (l 



1 - nr>, 1 



(72) 



where (p and <p are the energy obtained by the cell and 
transferred to the neighboring cell, respectively Hence, 
we can write the following kinetic equations for the en- 
ergy transfer: 



ri c A = n v {pi n c c - n%), 

ri c B = n v {pi n c c + n c A - n B ), 

n c c = n v (p\ n^ + n B - n c c ), 



(62) 
(63) 
(64) 



n° v = n v (pc n D - I) + pc ^V^n^n^) +rj(r,t) (65) 

Here, the discrete Laplacian A must be understood with 
the Neumann boundary conditions. The steady state cor- 
responds to the conditions, N c = 0. This leads us to the 
following driving conditions for the closed boundary: 



P C 3 M 

(P C 2+P C 3 )/P~ C , 
(p\+P2+P c 3 )/p~ c 



l/p c 



(66) 
(67) 
(68) 
(69) 



To perform the standard renormalization procedure 
described above and find the recursion relations, the gen- 
erating function method can be employed. To this end, 
we introduce the generating function for the relaxation of 
a cell on the closed boundary as the following polynomial: 



The critical exponent results to be r c = 1.239. The cor- 
rection to the bulk critical exponent due to the half-plane 
geometry was not presented before. The agreement of the 
bulk and open boundary exponents with exact values al- 
lows us to expect such an accuracy in the case of the 
closed boundary. 



VIII. DISCUSSION AND CONCLUSIONS 



In this paper we have presented the detailed appli- 
cation of the DDRG to the sandpile model. We have 
concentrated on the BTW model which we have studied 
using schemes of increasing complexity. In the simple 
scheme the sites are subdivided in three states (stable, 
critical and active) and the RG transformation acts on 
the energy transfer probabilities pi [ O ], The scheme 
is then extended in order to treat explicitly the four 
states probability densities n Q , which can be obtained 
sclf-consistently. The fixed point values of n a and pt are 
in good agreement with exact results. In addition, we 
compute the critical exponent r describing the avalanche 
size distribution. The result is in good agreement with 
numerical and analytical estimations and appears to be 
robust with respect to the different approximations. Fi- 
nally, we study the boundary scaling of the sandpile 
model, obtaining results in good agreement with exact 
results. The present analysis shows that the DDRG is 
a complete and systematic tool to study nonequilibrium 
systems in general and sandpile models in particular. 
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Po 


Pi 


P*2 


P*3 


p% 


p* =0.595 


0.091 


0.345 


0.379 


0.161 


0.024 



TABLE I. The fixed point probabilities for the energy 
transfer from the relaxing cell, including the probability po. 





Pi 


Pi 


Ps 


Pi 


Bulk 


0.295 


0.435 


0.229 


0.0414 


Open 


0.142 


0.417 


0.351 


0.0899 


Closed 


0.526 


0.394 


0.0799 





TABLE II. 
scheme. 



Relaxation probabilities in the extended RG 



»A 



nc 



Bulk 



RG 



0.0205 0.134 0.349 



0.496 



Exact [24 1 0.0736 0.174 0.306 0.446 



Open 
boundary 



RG 



0.0377 0.184 0.359 0.419 



Exact [25 1 



0.104 



0.217 0.316 



0.363 



Closed 
boundary 



RG 



0.0514 0.305 0.643 



Exact [25 1 



0.113 



0.318 0.568 



TABLE III. Height probabilities in the stationary state. 



RG 



1.248 



To 

1.486 



1.239 



Exact 



1.25 



1.5 



TABLE IV. Critical exponents r,r and r c for the 
avalanche size distribution in the bulk, open and closed 
boundaries, respectively. 
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FIG. 1. Example of the renormalization scheme for the relaxation dynamics. For details see the 
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FIG. 2. Samples of application of generating functions. 
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FIG. 3. Renormalizating the lattice we couple the two bonds connecting neighbouring blocks into the only bond of the new 
lattice. N u N 2 -> N, E 2 , E 3 -> S, S 3 , S 4 -> W, W 4 ,Wi S. 



1G 





FIG. 4. We show the four different types of the blocks for the cemichal reaction model and some relaxation schemes spanning 
them. The other schemes can be obtained from these figures by rotations. It is convenient for calculations to subdivide the 
relaxation processes in the block (d) into three parts shown in the dashed box. 
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